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I review recent results on QCD at high temperature on a lattice. Steady progress with staggered 
fermions and Wilson type fermions allow a quantitative description of hot QCD whose accuracy in 
many cases parallels that of zero temperature studies. Simulations with chiral quarks are coming 
of age,and togheter with theoretical developments trigger interesting developments in the analysis 



of the critical region. Issues related with the universality class of the chiral transition and the 
^ fate of the axial symmetry are discussed in the light of new numerical and analytical results. 

Transport coefficients and analysis of bottomonium spectra compare well with results of heavy 
I ion collisions at RHIC and LHC. Model field theories, lattice simulations and high temperature 

^ systematic expansions help building a coherent picture of the high temperature phase of QCD. 

The (strongly coupled) Quark Gluon Plasma is heavily investigated, and asserts its role as an 
CO inspiring theoretical laboratory. 
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1. From here to high temperature QCD 

High temperature QCD deals with phenomena related with the deconfinement transition to- 
wards Quark Gluon Plasma, and the Quark Gluon Plasma phase itself [1]. 

When heating up matter from room temperature towards temperatures which are best measured 
in MeV rather than in ordinary degrees we will first produce a gas of the lightest particles - a gas 
of pions. As temperature increases a more general gas of resonances, the Hadron Resonance Gas, 
is though to be a better approximation of the system. Beyond temperatures corresponding to the 
nucleosynthesis of the lightest nuclei we enter what we can consider a purely hadronic world. At 
a temperature of about 200 MeV - 2^ 210^^ K - hadrons cannot exist anymore, and a transition 
to the Quark Gluon Plasma takes place. The plasma survives till the Electroweak transition, at a 
temperature of O(IOO) GeV. All in all, the field of high temperature QCD spans about five orders of 
magnitude in energy, between the nucleosynthesis of light nuclei up to the Electroweak transition. 
As we - lattice theorists - know, many of the high temperature most significant manifestations are 
outside the reach of any analytic approach, making lattice studies mandatory. 

This paper describes the developments of this field after last year's Lattice review [2]. The 
material fits roughly in two main avenues : one more phenomenological line, which uses QCD, 
or approximation thereof, to describe mostly the physics of high temperature around the critical 
region, up to temperatures of about 600 MeV. High statistics, physical values of the quark masses, 
continuum limit are key ingredients here. One second line uses models or approximate calculations 
to learn about general aspects of the same phenomena. This second line reaches up to far higher 
temperatures, witnessing the deep theoretical interest of the yet largely unexplored physics of the 
Quark Gluon Plasma. 

Before closing this short introduction, one word about the experimental situation. During the 
first years of operation at RHIC there was a strong focus and motivation to explore the physics of 
small baryon densities - which, in the Gran Canonical approach we use, corresponds to the physics 
of ji/T < I., ji being the quark chemical potential. The current heavy ion collisions at the LHC 
create instead near zero density systems, ji/T « I while the search of the elusive endpoint of 
QCD and the study of other aspects of colder, and denser matter ji/T > \ will be investigated at 
RUN II at RHIC, and future experiments at FAIR and NICA. This review will concentrate then on 
the first aspect - the physics of zero baryon density. The current effort at non-zero baryon density, 
towards the understanding of the difficult, technical problems still hampering the simulations at 
large chemical potential, is discussed by Gert Aarts [3]. 

2. Thermodynamics : from hadronic Phase to QGP 

The transition from the hadronic phase to the Quark Gluon Plasma is a crossover for phys- 
ical values of the quark masses[2]. As such, the transition temperature is not uniquely defined 
and in past years a considerable effort has been devoted to accurate measurements of the pseu- 
docritical temperature identified by different observables. With staggered fermions there is now 
agreement among different estimates, when the same observables are being used. A key ingredient 
in resolving this issue has been a careful scale setting procedure. Still not completely converged 
before Lattice 2012 was the situation of the thermodynamics equation of state, namely the relation 
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Figure 1: Equation of state for A/j = 2 + 1 flavours, staggered quarks. The HotQCD coUaboration use 
HISQ/tree action, the Wuppertal-Budapest a stout action. The pion mass in either cases is close to its 
physical value. Note the different temperature ranges in the two diagrams (Ref. [6]). 

between pressure, energy and temperature. The phenomenological relevance is great, since the 
equation of state is an input to hydrodynamic calculations [4]. The quantity which is more easily 
accessible in a lattice calculation is the trace anomaly £ — 3P (often called interaction measure) 
which can be computed by differentiating the partition function [5]. It is customary to compare 
its behaviour to that of the Stefann-Boltzmann limit, corresponding to free fields: convergence to 
this limit is achieved only at very large temperatures, T > 10^^ GeV. This late convergence is taken 
as one of the evidences of the interacting nature of the Quark Gluon Plasma close to the critical 
temperature. 

2.1 Staggered fermions with physical values of u,d,s masses 

Only two collaborations have results with physical values of the up,down, strange mass, in 
either cases with staggered fermions. In literature and talks prior to Lattice 2012 there are reports 
of a discrepancy between HotQCD results on the interaction measure £ — 3P and the Wuppertal- 
Budapest stout continuum estimate. 

In fig. 1, from Ref. [6] we can see the plot of the trace anomaly (e — 3P)/T^ as a function of 
the physical temperature, as obtained from the two major staggered collaborations using stout [7] 
and HISQ staggered quarks [6] and physical quark up,down and strange masses. Once the scale 
is properly set by using /k, for a given value of Nt the results nicely agree with each other. The 
extrapolated data - according to the Wuppertal-Budapest collaboration - however lie significantly 
below the results for any finite Nt, suggesting residual discretization effects. As thermodynamics 
is now entering a precision era, perhaps it is useful to fix some significant value of the temperature 
and compare in details results - similarly to what we would do at T = 0. I have chosen T = 200 
MeV - in correspondence to the peak of the interaction measure - and T = 300 MeV, see Fig. 
2,left. Since for a fixed Nt the agreement is very good, I only show the finite Nt HISQ results, 
together with the stout extrapolated results. The relevance of the continuum extrapolation appears 
clearly in the plots - the extrapolation of the HISQ result depends crucially on which Nt one would 
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Figure 2: The interaction measure for two selected temperatures as a function of the l/N^ : results at fixed 
Nt compare favourably, some residual issues with the extrapolation procedure are being settled. On the 
right, the overall status of the continuum results for T = 300 MeV (Numerical results communicated by the 
Authors). 



include. In this sense, the resolution of the residual small discrepancy is narrowed down to the 
issue of Nt extrapolation, and either collaborations are steadily improving their results. The same 
Fig. 2,right, shows the dependence of bulk thermodynamics on the matter content, to be discussed 
below. 

2.2 Wilson fermions thermodynamics 

New results for the thermodynamics with Wilson fermions have been presented by the WHOT- 
QCD collaboration[8] and by the tmft collaboration [9]. In either cases the pion mass is larger than 
400 MeV. 

The WHOT-QCD collaboration presented results for 2 + 1 flavors, a physical value of the 
strange quark mass, and a largish pion mass. Some results are shown in Fig. 3. The Wilson action 
was non-perturbatively improved, and the gluonic part used Iwasaki improved glue. Results for two 
flavors of twisted mass Wilson fermions at maximal twist are available as well, with a pion masses 
of about 400 MeV, and larger, see Fig.4. The results used the tree-level correction advocated by 
the Wuppertal Budapest collaboration, and the extrapolation to the continuum limit appears to be 
reasonably under control. 

2.3 Effects of the charm on the EoS 

One step ahead in the matter content was taken by the MILC collaboration and by the Wuppertal- 
Budapest Collaboration who presented results for A// = 2 + 1 + 1. Contrary to naive expectations, 
and in agreement with the predictions of effective theories[l 1], there is an early onset of the charm 
contribution at about T = 350M^y. The MILC collaboration [10] determined the QCD EoS with 
2+1 + 1 flavors using the HISQ action. Their study was then carried out on lines of constant 
physics with mi/ Ms = 1/5, with physical strange and charm quark mass. Above Tc the precise 
value of the light quark mass should not be important. The preliminary data have been obtained 
with Nt=6, 8, 10 and 12, without any attempt at a continuum extrapolation yet. Fig. 5, left shows 
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Figure 3: Equation of state for A^/ = 2 + 1 improved Figure 4: Equation of state for Nf = 2 Wilson fermions 
Wilson fermions. mj^/mp = 0.63 (Ref. [8]). at maximal twist, 400 MeV (Ref. [9]). 
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Figure 5: Charm effects on the Equation of State: MILC results for the trace anomaly for A^j = 2 + 1 + 1 
using HISQ quarks along lines of constant physics(left) and comparison of results with and without dynam- 
ical charm(right). This study was carried out with physical charm and strange mass. The comparison of the 
left and right plots allows the identification of the effect of the charm mass on the EoS (Ref. [10]). 



the MILC preliminary result for the trace anomaly e — 3P , and in Fig. 5, right, the same is shown 
but with either without the valence charm quark contribution. 

The Wuppertal-Budapest collaboration[12] presented further results on charm effects on the 
EoS. Their results were also successfully compared with High Temperature perturbation theory, 
down to r = 350 MeV, and with a Hadron Resonance Gas model, which offers a good description 
of the system up to temperature of about 150 MeV. The outcome of this comparison is presented 
in Fig. 6. As discussed at length in previous studies, any simple analysis is expected to break down 
in the vicinity of the transition, where the critical or pseudocritical behaviour cannot be captured 
by an analytic equation of state. Fig. 6 thus clearly identifies the interesting region where a lattice 
analysis is mandatory. An overview of the mass sensitivity for T = 300 MeV can be seen in Fig. 2, 
right, where all available continuum results are shown as function of the pseudoscalar mass. 



3. The transition region - symmetry aspects 

In the limit of massless quarks the QCD Lagrangian enjoys a global chiral symmetry SU (A^/)l ^ 
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SU {Nf)R (7a(1) ® CA^(l)- The vector C/y (1) transformation corresponds to baryon number and 
the axial C/a(1) which would promote the symmetry to U (Nf) ® U (Nf) is a symmetry of the clas- 
sical theory but not of the quantum theory. The symmetry breaks down to SU {Nf) ® 6V(1) at low 
temperatures. Since up and down quarks, and to some extent the strange quarks, are not too heavy, 
relics of the critical behaviour can still be observable even in the crossover region for physical 
values of the quark masses - the connection being made within the framework of the magnetic 
equation of state. 

It is speculated that the breaking of the 
Ua{1) symmetry is stronger at zero tempera- 
ture, or, equivalently, that C/a(1) symmetry 
is effectively restored at high temperature. 
The fate of the C/a(1) effective restoration 
is closely related with the issue of the order 
of the chiral transition with two massless fla- 
vors, leaving open the possibility of either a 
first or second order transition depending on 
the intensity of the residual C/a(1) breaking 
atrc[15, 16]. 

In view of the relevance of the symme- 
tries, it is obviously extremely desirable to 
have fermionic formulations which respect 
chiral symmetry. Remember now that stag- 
gered and Wilson fermions break the C/a(1) 
symmetry, and that chiral symmetry is also 
broken in the Wilson formulation, while an (7(1) symmetry (unrelated with the C/a(1) symmetry) 
is preserved by staggered fermions, which however do not allow a local definition of (7a(1) sym- 
metry. The computationally expensive overlap quarks have an exact chiral symmetry. Domain wall 
fermions approximate it with a precision which is in principle arbitrary, the residual chiral sym- 
metry breaking being characterised by an additive renormalisation to the quark mass rrires which 
vanishes in the limit of infinite 5-th dimension: rriresiLs) = ci ^-r-^ + C27-. 

New results with Wilson fermions with 2+1 and 2 flavours, as well as results with chiral 
fermions were presented at this conference, with emphasis on two interrelated open issues: the fate 
of the C/a(1) symmetry at the phase transition, and the universality class for Nf = 2. 

3.1 C/a(1) 

The physics of the anomaly is particularly intriguing, given its role in the universality class of 
the transition [15, 16]. Moreover, should the anomaly not be restored at the transition, one would 
have the opportunity to disentangle its effects from those of chiral symmetry. 

Patterns of symmetries are reflected in the interplay of scalars susceptibilities, according to 
the well known MILC plot [13, 14]. The order parameter for 5C/(2) 5C/(2) is related to the 
longitudinal susceptibility <\j/yf>/m = X7i^ ^ candidate order parameter f or Ua(1) is Xti — Xo- 
It is however true in general that if an observable <^ is not invariant under a given symmetry 
< >= is a necessary, but not sufficient condition for y to be exact. This means that the 
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Figure 7: Direct observation of the degeneracy in the scalar sector, right above Tc, left. The eigenvalue 
distribution for several temperatures, hinting at the restoration of the Ua{1) symmetry. Overlap fermions, 
two flavors (Ref.[22]). 

candidate order parameter for the C/a(1) symmetry Xn — Xo might well be zero even if (7^(1) is still 
broken. 

Chiral observables are related to the low-lying eigenmodes of Dirac operator: 

/2m 

Xn-Xs = J dXp{X)j-^^^ (3.2) 

so the fate of the C/a(1) synmietry can be read off the shape of the eigenvalue density p(A). The 
various possibilities are extensively discussed discussed in a recent paper by the HotQCD Collab- 
oration [21]. 

A new theoretical analysis has been presented by Aoki[17, 18] and collaborators, who find 

lim < p^(A) >m= lim < >^ — + 0(A4) (3.3) 

i.e. no constant contribution to the spectrum, which implies at high temperature 

lim(Z;r-ZcT) = 0. (3.4) 

when < xj/Xj/ >= 0. In other words, for any Nf, two point correlators - hence susceptibilties - 
would be blind to the breaking of Ua{1) when chiral symmetry is restored. Note that older results 
(refs. [19, 20]) reached a similar conclusion on the non-observability of C/a(1) non restoration in 
two point functions, but only for Nf>3. The new analysis [17] requires a careful consideration of 
the order of the limits and suggests that early work was inaccurate. 

Results presented by JLQCD[22], with two flavors of overlap fermions and physical values of 
the pion mass, apparently realise the scenario proposed by Aoki[17] indicating the degeneracy of 
all the scalar correlators in the chiral limit above T^, which is of course a stronger statement than the 
degeneracy of the susceptibilities, see Fig. 7, left. Moreover, a gap in the eigenvalues' distribution 
opens up right above Tc, as shown in Fig. 7, right. The conclusion still needs to be corroborated by 
a better control of systematics. 
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Figure 8: Domain wall fermions via-a-vis staggered results (top left). The scale independence for domain 
wall fermions (bottom left). The eigenvalue spectrum (right) favours a non restored Ua{1) symmetry be- 
yond the QCD transition, however the largish masses and the residual breaking leave room for a different 
conclusion (Ref. [23]). 



The same issue was investigated with domain wall fermions. New results with 2+1 flavours 
were presented at this meeting[23] on the 32^ x 8 lattices by RBC/LLNL, and ensembles with 
physical quark masses on the 32^ x 8 and 64^ x 8 are being planned from HotQCD. The behaviour 
of chiral susceptibilities which can be seen in Fig.8, left top, is on the overall convincing, with small 
deviations from the existing staggered results. Fig.8,left bottom, shows the ratio of dimensionless 
ratios at j8 = 1.75, demonstrating violation of the order of only 5 % from an ideal scaling. The 
spectral densities are shown in Fig.8, right, for a temperature above, but close to T^. The distribution 
shows a clear peak near zero ("near zero modes") and a linear behaviour extending down to zero. 
This may correspond in the continuum limit to p(A) = ci5(A) + C2A. at a variance with Aoki's 
scenario, and implying a residual C/a(1) breaking with Xti — X5 7^ ^ the chiral limit. 

All these analysis require an extrapolation to the chiral limit, hence studies with different 
masses - which are in progress - are needed to draw definite conclusions. However a detailed study 
with HISQ fermions[24], 2+1 flavors, pion masses ranging from 115 MeV till 230 MeV leaves 
the conclusion open, with a careful discussion of the systematic effects affecting a proper analysis 
of the eigenvalue shape.The results for the distributions are presented in Fig. 9. Even at the highest 
temperatures there is no evidence of a gap around zero. More than any comment the shape of the 
eigenvalues clearly demonstrate the uncertainties affecting the conclusions. 

3.2 Further studies around the critical temperature 

Further studies with chiral fermions were presented by the TWQCD Collaboration , using two 
flavours of Domain Wall fermions [25] and by the Wuppertal-Budapest collaboration[12, 26], who 
use two flavours of overlap fermions, with topology fixing a la JLQCD. These studies should be re- 
garded as still preparatory since the pion masses exceed 300 MeV, and the analysis concentrates on 
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Figure 9: Eigenvalue distribution from HISQ fermions: the shape determines the fate of ^7^(1) symmetry 
which remains open. The left diagram shows temperature dependence of the Dirac spectrum at temperatures 
in the range 173.0 — 331.6 MeV > Tc. The right diagram shows the volume dependence of spectral density 
at the highest temperature. The volume dependence is quite small and there remains a tail even at this 
temperature (Ref. [24]). 
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Figure 10: The crossover region with overlap fermions, A// = 2, ruji = 350 MeV (Ref. [26]). 




Figure 11: The renormalised chiral condensate, renormalised Polyakov loop and strange quark susceptibility 
from the Wuppertal-Budapest simulations with 2+1 flavors of Wilson fermions (Ref. [27]). 



basic observables. The Wuppertal-Budapest collaboration presented -see Fig. 10 - a comparison 
with continuum staggered results which is at first sight satisfactory[12, 26]. The transition region 
was also studied with Wilson fermions, 2+1 flavors, including continuum extrapolation [27, 29]. 
The pion masses range from 125 to 275 IMeV, and studies with smaller masses are in progress. 
Again the focus is on the agreement with staggered quarks, see Fig. 11, which so far looks very 
nice, within the largish errors. New results for two flavors of non-perturbatively 0{a) improved 
Wilson fermions have been presented [28] working along lines of constant physics with ruji 290 
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Figure 12: Vector Correlators in the plasma, from simulations with two flavors of improved Wilson fermions 
on large lattices (left), and summary of results (right)(Ref. [28]). 



MeV. The picture, Fig. 12, shows the axial vector correlator normalised to the free vector correla- 
tors. The right-side picture shows the collection of screening masses themselves - the degeneracy 
in the vector sector which anticipates the one in the scalar sector might suggest a persistence of 
Ua{1) breaking, although it could also be that the near degeneracy in the vector sector comes from 
the largish masses' values. 

3.3 The universality class for Nf = 2 

As noted above, the order of the transition for two flavor QCD is a dynamical issue which 
depends on the magnitude of the residual breaking of Ua{1) at the transition point. 

The situation can be depicted as in Fig. 13, left: the order might well change along the 
transition line even if the synmietries are the same on both sides. Since quarks chemical potentials 
do not alter the chiral symmetries in the Lagrangian one can study a possible realisation of this 
scenario by varying /i, away from the condensation phases, and see whether one can observe first 
and second order transitions. This strategy is pursued [30, 31, 32] with imaginary chemical 
potential. Fig. 14 shows the results with imaginary baryonic jU[30], and Fig.l5 shows the results 



The order migh change 




rriTrc 



Figure 13: Sketchy view of a possible behaviour of the transition line: the order might change along the 
line even when the symmetries of the two phases are the same(left); the region of the first order transition is 
characterised by the endpoint at finite mass, in the universality class of the Z2 Ising model (right). 
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Figure 14: Search for universality class of the Nf = 2 theory : the position of the endpoint of the first order 
transition for imaginary /i versus the bare mass (the exponent in the right diagram is motivated by the Z2 
universality class associated to the endpoint). The transition becomes of second order when the endpoint 
hits zero, and based on these plots the authors suggest that this might happen for a positive jl^. The results 
have been obtained for Nt = 4 unimproved staggered fermions hence they are far away from the continuum 
(Ref. [30]). 



for imaginary isospin jU[31]. The authors argue that following the fate of the endpoint of the first 
order transition, sketched in Fig. 13, right, is a practical way to decide the order of the transition at 
jU = 0. In particular one might want to locate the turning point between a first and a second order 
transition: since first order transitions were observed for a non-zero imaginary chemical potential 
for A// = 2, a turning point for real ji implies that the phase transition at zero chemical potential is 
of first order [30, 31]. Current extrapolations suggest that this is realised for Nt = 4, for unimproved 
staggered fermions, hence far from the continuum limit. 

The transition for two flavor close to the continuum limit is being investigated with Wilson 
fermions [9], and some results are shown in Fig. 16. The plot depicts the current situation with 
Wilson fermions, maximal twist and improved gauge sector. The extrapolation to the chiral limit 
can be carried out either with a second order ansatz, or assuming a critical endpoint smaller than 
300 IVleV. The data do not discriminate between these two possibilities. This implies that the order 
of the Nf = 2 transition in the massless continuum limit is still undecided, and suggests that the 
endpoint of the first order transition at /i = 0, if any, is for ruji < 300 IVleV. 
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Figure 15: Direct observation of a first order transition at non-zero imaginary isospin density. If this signal 
would persist till zero chemical potential the Nf = 2 transition would be of first order (Ref. [31]). 
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for Nf = 2, with chiral extrapolation as moti- 
vated by the plausible universaHty classes (Ref. Figure 17: The phase diagram in an enlarged 
[9]). parameter space (Ref. [34]). 



4. The transition region - confinement, gauge dynamics, external fields, G term 

The discussions in the previous Section focused on the fermionic sector of the theory, and the 
associated chiral symmetries. The "obligatory" step at this point is to remind ourselves that gauge 
dynamics and confinement are crucial, yet mysterious, aspects of the transition, in particular due 
the lack of a bona-fide order parameter for deconfinement. 

Abelian gauge theories are an important exception : in these theories the confinement transition 
is understood as a condensation of topological defects - in the case of a four dimensional i7(l), for 
instance, strings of magnetic monopoles [33]. 

Ogilvie proposes to search for a continuous path between non-abelian theories and abelian 
ones by suitably enlarging the parameter space, [34]. The considered phase diagram is depicted in 
Fig. 17. By supplementing the Yang-Mills SU (2) action with a double trace term, one can show 
that the confined phase of pure SU (2) and the high temperature confined phase of the extended 
model are continuously connected; the latter is amenable to a semiclassical treatment and in turn is 
continuously connected to the confined phase of the abelian model. 

A more direct interpretation of the transition in terms of magnetic monopoles has been recently 
proposed by Shuryak [35]. In Fig. 18 lattice data for different Nf are used to estimate the strength 
of the coupling at the transition, which is then compared with the coupling at the IRFP[36, 37]. 
In particular the data marked with a diamond were obtained with the same action,and different 
number of flavors, and can be meaningfully considered together [38]. The coupling at the tran- 
sition becomes stronger and stronger with increasing Nf: the proposed explanation is that light 
fermions can occupy the (chromo)magnetic monopoles, making them unsuitable for Bose Einstein 
condensation. Lattice simulations measuring magnetic monopole properties can confirm this pro- 
posal, which,if true, would lend support to confinement scenario based on monopole condensation 
in QCD. A further related observation is the emergence of the conformal window of QCD for large 
number of flavours - reviewed by Joel Giedt[39] - as the zero temperature limit of a super-strongly 
interactive Quark Gluon Plasma. 

It has been proposed that fluctuation of topological charge might play a role at the transi- 
tion, producing measurable P and CP violations in heavy ion collisions. In the QCD Lagrangian 
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Figure 18: The transition as a function of matter content (left) and the associated coupling according 
Results from various lattice groups (Ref. [35]). 
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Figure 19: News on the 6 vacuum: the continuum extrapolation of the slope of the critical temperature from 
lattice simulations at imaginary 0(left, Ref. [40]), and from model studies (right, Ref. [42]). 



P and CP violating terms enter via the 9 term OQ{x), Q{x) being the topological charge density. 
Despite the experimental bound < 10~^^ it is still conceivable that such term affects the criti- 
cal temperature. New studies on the 0-dependence of the deconfinement temperature have been 
presented[40, 41, 42], the former numerical and the latter based on a model analysis. In both cases 
T{9) was parametrised as T{9) = Tc—RqO^, and the coefficient of the quadratic coefficient /^^ has 
been estimated. The results are shown in Fig. 19. Both studies agree on the sign of Rq term, as well 
as on its very little magnitude, with a substantial agreement within the largish systematic errors. 

Finally, an interesting line of research is the effect of an external field on the phase transition[43]. 
As magnetic fields as large as 10^"^ Tesla can be produced at the QCD transition a lattice effort has 
started to compute the effects of strong magnetic fields on the QCD transition, and on the magnetic 
catalysis - a possible enhancement of chiral symmetry due to the effects of the external fields. Al- 
though intuitively an external magnetic fields should favour chiral breaking, model analysis give 
non-univocal results, and, again, lattice calculations are mandatory to settle the issue. Details of the 
dynamics are still under discussion, particularly close to the chiral limit and I refer the interested 
reader to a recent review for details [44]. Recent results have been obtained for two [45] and three 
colours [46]. 
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5. Fluctuations 

Lattice QCD simulations are being carried out in a GranCanonical ensemble. This allows fluc- 
tuations of conserved quantities - baryon number, electric charge, strangeness - for instance, the 
baryon number should be zero on the average, but it can fluctuate configurations by configurations. 
In ordinary conditions fluctuations are quite small: for instance one can compute the distribution of 
baryon number in ChPT and find that it is described by a narrow Gaussian [47, 48] while when the 
temperature increases the distribution broadens and might change shape[47, 48]. The fluctuations 
thus reflect properties of the thermal medium and of the critical point[49]. 

In general, a synthetic way to describe a distribution is via its momenta 

/ + 00 
X'~f{X)dX (5.1) 
-oo 

The momentum generating function M((^) =< exp{^X) >= < > is a simple way to 

combine momenta into a single expression and we recognise that the momenta are linked with the 
coefficients of the Taylor expansion of the generating functional < >= [^M((^)]^^o- Consider 
then the pressure p = T /Vln^{T^ii) and its Taylor expansion in powers of ji/T around /i = 0: 

^ = E.,(n&» (5.2) 

and we recognise that the zero density fluctuations of the quark number density — ^^^fl^^ 
given by the coefficients Cn{T). Similar reasoning holds for the other conserved charges. Results 
for the momenta have been successfully compared with Hadron Resonance Gas till temperatures 
very close to Tc [50, 51]. 

At the time of the meeting new analytic results from A. Vuorinen and collaborators [52] be- 
came available, followed by a publication [53]. The excellent agreement between analytic and 
numerical results - which can be appreciated in Fig.20 - is in part due to the elimination of the 
pure gauge component (by taking the derivatives w.r.t to /i.) Note also in the same figure the very 
good agreement between the results of the Wuppertal-Budapest collaboration and HotQCD: these 
measurements have been performed on the same set of gauge fields discussed above and the agree- 
ment here again corroborate the optimistic view that the residual minor disagreements are almost 
solved. 

A completely different approach to the dynamics of fluctuations is offered by the flux tube 
model presented by Patel[54]: the idea is that the fluctuations should be 'visible' in the gauge 
fields, and one can conceive, within the flux tube models, one way to visualise the gauge field 
dynamics in position space. The 'rules of the game' in a first approximation are simple: quark 
and anti-quark are bond by a string in a QQ state. The string can break via the formation of a 
light-heavy meson, QQ qQ-\-qQ , or in the plasma as QQ Qq + qqQ, see Fig 21. In this 
way one might hope to study multi-correlation particles which are useful for phenomenology, and 
highlight the microscopic structure of the medium, for instance by distinguishing a dilute liquid 
from a cristal-like structure. 
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Nf=2+1, perturbation theory vs. lattice data; Andersen, Mogliacci, Su, Vuorinen (in preparation). 
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Figure 20: Lattice results for the baryonic susceptibility vis-a-vis analytic studies. Note the agreement 
between different lattice calculations, and with the 4 loops dimensionally reduced theory (Ref. [52]). 



6. Quark Gluon Plasma - a new world 

The physics of the Quark Gluon Plasma is of course addressed by most of the studies described 
before. Here -somewhat arbitrarily - 1 collect those studies which are motivated by the desire of a 
deeper and general theoretical understanding of the physics of this unusual state, without reference 
to specific models. 

One popular way to learn about general features is to enlarge the parameter space of interest: 
a large number of colours, or a large number of flavors. For a large number of colours Nc - 
reviewed by M. Panero[55] - we have one very interesting finite temperature result[56]. Panero 
and collaborators have studied the behaviour of the Polyakov loop in the deconfined phase of 
a SU {Nc) Yang-Mills theory. Different representations and different Nc exhibit a near universal 
behaviour once geometric scaling factors are taken into account, see Fig. 22. This observation - 




Figure 21: Describing fluctuations: the flux tube model by Patel (left, Ref. [54]). The Polyakov loop 
analysis by Gattringer (right, Ref. [59]). 
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SU(6), N=5 



SU(4), fundamental representation 



fundamental 




Figure 22: The left diagram shows the bare loop for SU (6), and different representations, rescaled assuming 
Casimir scaling-the universality of the results support the assumption. The right-hand diagram shows the 
same renormalised loop in the fundamental representation for SU (4) -in both cases the behaviour is very 
similar to the one observed in SU (3) (Ref. [55]). 

so far restricted to pure gauge systems - lends support to large Nc approaches to the analysis of the 
high temperature phase of QCD. 

It is believed that as chiral symmetry is no longer spontaneously broken, the only remaining 
scale in the massless model is the temperature. Would the eigenvalue spectrum of QCD give 
information on any remaining scale? Kovacs[57] has shown that this is indeed the case - the 
equivalent of the Thouless energy was found to scale as the temperature itself. 

Attempts at modelling the Quark Gluon Plasma by use of Polyakov loop dynamics has pro- 
duced the observation that Z3 metastable states might become important at T > 750 MeV [58]. 
Such metastable states with dynamical fermions might have an influence on bubble formation and 
on the system's dynamics in general. 

How to detect directly the dynamics of the Polyakov loop? Diakonov, Gattringer and Schadler 
[59] propose a free energy for the untraced Polyakov loop X 



where X is defined by 



and 



I Nt 

— Y UAx,t)^TrX 



< 



5>r,x,^'=-^lnZ(r,X,i80 



(6.1) 



(6.2) 



(6.3) 



This definition, so far tested in the pure Yang Mills case, see Fig. 21, right, compares very well 
with high T perturbation theory and opens the way to a rigorous approach to the study of Polyakov 
loop dynamics. 

7. Quark Gluon Plasma confronts experiments 



How far we can go with comparing our results with experiments? For orientation, we know 
that during Au-Au collisions at 200 GeV at RHIC the system reached up temperatures in the range 
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340-380 MeV, the current LHC runs - 2.76 TeV - explored up to 420-480 MeV with possibly hot 
spots of about 500-600 MeV. LHC at full luminosity will reach about IGeV - six times Tc [60]. 
The recent CMS results [61] are a good starting point for the discussion: the invariant mass plot 
measured in pp shows the expected three bumps in correspondence to the fundamental T state and 
the first two radial excitations. When the same distribution is being measured by colliding lead 
ions, the lowest peak is almost stable, while the two higher peaks have almost disappeared. A 
natural thing to do is to compare these observations with the heavy quark spectrum computed on 
the lattice - a notoriously difficult task, see eg. Ref [62]for a recent comprehensive overview. 

In early times a lattice computed potential was used as input for the Schroedinger equation for 
heavy mesons. Later on, issues were risen concerning the appropriate potential to use, and most 
important it was finally realised that the correct potential should include an imaginary component 
to account for the broadening of the state [63, 64]. Alternative calculations of the potential using 
the wave function as an input to the Schroedinger equation, and a reverse engineering to extract the 
potential itself have been presented at this meeting as well[66], and the results can be seen in Fig. 
23. 

At the other extreme - w.r.t. quantum mechanical, non-relativistic studies - we have field 
theoretic, relativistic calculations of the spectral functions:the euclidean charmonium correlator 
7 is computed[67] - in this case on quenched configurations -, and the spectral function is then 
reconstructed by use of by now standard Maximum Entropy Method analysis (MEM) [68]. 
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Figure 23: A new strategy for the potential analysis: Interquark potential obtained by reverse engineering 
the Schroedinger equation after calculations of the wavefunction (Ref. [66]). 




Figure 24: Charmonium spectral functions from a full relativistic calculation (Ref. [67]). 
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Figure 25: Momentum dependence of spectral functions in the plasma, as seen in the mass and width. 
Position of the ground state peak M{py)/M{0) (left) and the upper limit on the width of the ground state 
peak, normalized by the temperature, T/T (right), as a function of the velocity squared (v^) in the vector (T) 
channel (Ref. [70]). 

All these approaches give results broadly consistent with each other, and with intuition. How- 
ever a solid understanding of the systematics is still lacking. 

An improvement to the implementation of the Maximum Entropy Method has been presented 
by A. Rothkopf [65]. The method proposes to use an extended search space, and preliminary 
successful results have been obtained by use of mock data: the numerical behaviour is more stable, 
and the resolution is satisfactory. 

One further approach consists of evaluating NRQCD propagators over dynamical gauge fields: 
the reconstruction of the spectral functions is made easier[69] as in this case the spectral relation 
reduces to 

G(T)= r ^exp(-(o'T)p(a)') , (7.1) 

J-2M 71 

As a result, the spectral function is merely the inverse Laplace transform of the correlator, and 
all problems associated with thermal boundary conditions are absent. Remarkably, the NRQCD 
bottomonium spectral functions [69] show a suppression of excited states consistent with the CMS 
results [61]. The new NRQCD results [70] presented this year include a systematic study of mo- 
mentum dependence, as well as a study of the mass dependence[71]. Successful comparison with 
the predictions of effective models has started, together with a more detailed analysis of the sys- 
tematics. 

Directly related with the computation of spectral functions are the transport coefficients, see 
e.g Ref. [72] for a comprehensive review. New results appeared this year concerned heavy quark 
diffusion. In the framework of a Gaussian and Markovian Langevin dynamics, the diffusion 
constant D describes the mean square displacement for times t larger that the relaxation time 
< |A(x)p >= 6Dt. In the strong coupling limit of the SUSY Yang-Mills plasma one finds, us- 
ing holographic duality 
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Figure 26: Heavy quarks diffusion coefficients : impressive agreement with experimental results, corrobo- 
rating a strongly interactive quark gluon plasma - the limiting case is shown as a short segment on the right 
(Ref. [76]). 



Gubser[73] suggests A 2^ 6;r in QCD, concluding D 2^ l/(27rr). In a relativistic setup, transport 
coefficients are defined as 

p(co) 



lim 



(7.3) 



0)^0 CO^ 

where p{(o) is a spectral function in the appropriate channel[74]. For instance D has been extracted 
from the charmonium spectral functions described above[67], with results remarkably close to the 
one obtained in strongly coupled SUSY, eq. 7.2. 

The diffusion coefficient D can also be extracted from pure gauge correlators[75]: lattice 
results for this have been presented by Datta[76]. Again the result is close to the one of eq.7.2, 
hence in remarkable agreement with those from Ref. [67]. 

The comparison with experimental results is shown in Fig.26. Needless to say, it will be of a 
great interest to see the new results from LHC, at higher energies, when the plasma is expected to 
be less bound, and the transport coefficient should correspondingly increase. 

To appreciate the quality of the lattice results one should keep in mind that the LOPT results, 
see again Fig. 26 - to be achieved in the extremely high temperature limit - would be 

671 



D 



LOPT 



^80 



(7.4) 



g'^Tlog{2T/mD) 

On the scale of the huge range of allowed D values, from the extreme strongly coupled to the free 
plasma, the agreement between the two lattice calculations, carried out with different strategies and 
actions, is impressive. 



8. Summary 

A precision era for high temperature QCD has started, with results whose accuracy is often 
comparable with the one achieved at zero temperature. Staggered bulk thermodynamics in the 
continuum and at the physical point, well into the current LHC region is in good control, and 
residual discrepancies seem minor. Wilson fermions are becoming competitive, with results in the 
continuum limit (albeit still with largish masses) for bulk thermodynamics, and in the critical re- 
gion. New analytic studies in the high temperature phase compare remarkably well with numerical 
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results, down to unexpectedly low temperatures. Chiral fermions are coming of age, with new, 
continuously developing high quality results for domain wall and overlap fermions. Good chiral 
properties together with new theoretical insight trigger further activity on chiral and axial symme- 
tries and their spectral properties. Details of QCD dynamics prove to be important and delicate 
close to the phase transition: the fate of the axial anomaly, the order of the transition for 7V=2,the 
response to a magnetic field, and the role of a therm are still unclear. New lattice results for 
transport and quarkonia have been presented which compare well with RHIC and LHC experi- 
ments. The (strongly coupled) Quark Gluon Plasma is heavily investigated, either in QCD and by 
use of deformed QCD and other model field theories, and confirms its role of inspiring theoretical 
laboratory. 
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